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Discovery of an eclipsing X-ray binary with a 32.69 hour period in 
MlOl: an analog of Her X-1 or LMC X-4? 
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Kip Kuntz^ 

ABSTRACT 

We report the discovery of an eclipsing X-ray binary in MlOl, the first such 
system to be discovered outside the Local Group. Based on a sequence of 25 
Chandra observations that sample a wide range of orbital phases, we find a period 
of 32.688 ± 0.002 hours, which we interpret as an orbital period. The folded light 
curve exhibits an eclipse lasting about 8 hours, suggesting a compact orbit in an 
nearly edge-on configuration. The X-ray binary has an average luminosity of 
Lx(0.3 — 8keV) 1.3 x 10^® erg/sec, with only one out of the 25 observations 
significantly lower in flux than the average light curve. The presence of the eclipse 
and the ~ 1.4-day orbital period suggests that this source is an analog of the 
well studied eclipsing X-ray binary pulsars Her X-1 or LMC X-4. Combining 
the Chandra data and the HST ACS/WFC images, we have identified several 
possible optical counterparts, including an 05-03 star with V = 25.0. Follow-up 
optical monitoring observations should be able to identify the donor and further 
constrain the orbital properties. 

Subject headings: Galaxy: individual(MlOl) — X-rays: eclipsing binary 

1. INTRODUCTION 

The detections of X-ray binaries have recently begun in galaxies outside the Local Group 
with the advent of the Chandra X-ray observatory. Owing to its unprecedented sensitivity 
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and spatial resolution, it is possible to resolve individual X-ray point sources in these galaxies. 
In the past six years, Chandra has studied about one third (~70/220) of the galaxies with 
1^25 > 1' within 10 Mpc (Liu 2005); these galaxies span the whole spectrum of galaxy types, 
including spirals (e.g., M81, Tennant et al. 2001; M51, Terashima & Wilson 2004), ellipticals 
(e.g., NGC1399, Angehni et al. 2001; NGC4697, Sarazin et al. 2000), starburst galaxies (e.g.. 
Antenna and M82, Zezas et al. 2000), and colhding galaxies (NGC7714/7715, Smith et al. 
2005). The distributions and luminosity functions of the X-ray binaries in these galaxies, as 
compared to those in the Milky Way (Grimm et al. 2002) and M31 (Kong et al. 2003), show 
significant dependence on the galaxy type and star formation properties. 

An ultra-deep Chandra ACIS survey (PI: Kuntz) has been executed for the face-on spiral 
MlOl (NGC5457) at a distance of 6.8 Mpc (Freedman et al. 2001). Twenty four Chandra 
observations were tiled to achieve a full coverage of MlOl, with a cumulative exposure of 
one mega second within 4' of the nucleus. Such a deep exposure can reach a detection 
limit of 4x10^^ erg/sec, only a factor of two more luminous than current Chandra studies 
of M31 (Kong et al. 2002). This survey allows detection of most low-mass X-ray binaries, 
most high-mass X-ray binaries, most young supernova remnants and many middle-aged ones, 
super-bubbles and giant Hll regions. The observations were carried out over the year from 
January 2004 to January 2005, allowing studies of the variability and transient behaviors of 
these sources. 

Here we report the discovery of a 32.688±0.002 hour period for an eclipsing X-ray binary 
from the MlOl ultra-deep survey and an early observation in 2000. The Chandra observations 
and the search for periodicities are described in section 2. Additional optical observations 
with HST ACS are described in Section 3. In Section 4, we describe the constraints the 
32.69 hour period places on the accretor and donor masses, and discuss the prospects to 
monitor and identify the donor at optical wavelengths. This X-ray source has been detected 
previously in ROSAT HRI and PSPC observations and named as M101-X7 (Liu & Bregman 
2005; H36/P21 in Wang, Immer & Pietsch 1999), but the periodicity could not have been 
detected given the low photon counts and the small number of observations. 

2. Chandra Observations 

M101-X7 was observed in 25 Chandra observations spanning five years (Table 1), in- 
cluding 24 observations from the MlOl ultra-deep survey and an earlier observation (ObsID 
934, 2000-03-26). All the observations were downloaded from the Chandra Data Archive, 
and processed with CIAO 3.3.0.1 and CALDB 3.2.1. WAVDETECT was run with the scale 
set of (0'.'5, I'.'O, 2'.'0, 4'.'0) to detect point sources on the X-ray images. 
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The source was detected in 25 observations at an average location of R.A. =14:03:36.0, 
Decl.=54:19:25, thus named as CXO J140336.0+541925 following the CXO naming conven- 
tion. The source was at large off-axis angles in most observations, and was elongated at 
different elongation angles due to the different positioning of the aim points. The merged 
profile from all observations was symmetric and could be fitted by a Gaussian with FWHM 
of 1'.'7, consistent with the point source diameter (1'.'7) for enclosed energy fraction of 50% 
at an off- axis angle of 5 arcminutes. 

The light curve was constructed to demonstrate variations during individual observa- 
tions. The photons (0.3 kcV <E< 8 kcV) were extracted from the source region enclosing 
95% of the total photons reported by wavdetect. The time bin was adjusted to have an aver- 
age of >15 counts per bin and no more than 128 bins. Among the 25 light curves, emphases 
were put on a group of 18 "good" hght curves for the 18 observations with more than 100 
photons, with 17 observed in 2004 during the MlOl ultra-deep survey, and one (ObsID 934) 
observed in early 2000. Among the 18 "good" light curves, 12 light curves showed apparent 
variations suggestive of periodicities, with apparent minima in seven of them (Table 1). 

Two techniques were apphed to search for possible periods. The first was the Lomb- 
Scargle method devised for unevenly spaced data (Lomb 1976; Scargle 1982). The Lomb 
periodogram was computed for the bins in the 18 "good" light curves (Figure 1), which 
showed significant peaks at 32.69 hours, 22.44 hours, 41.75 hours, and 16.10 hours. The 
probabilities for these peaks to originate from random fluctuations of photons in the absence 
of true signals arc less than 10"^'^, suggestive of the presence of true signals. This, however, 
docs not mean all peaks in the Lomb periodogram are true periods. Indeed, some peaks 
may be just period aliases. Another method applied to the data was the phase dispersion 
minimization method as demonstrated in Stellingwerf (1978). A rough search of periods in 
the range of 5 hours to 50 hours with steps of 0.045 hours revealed significant phase dispersion 
(6) minima at 32.72 hours, 44.87 hours, 33.53 hours, 31.91 hours, and 22.82 hours (Figure 
2). Searches were then carried out with finer steps of 2 x 10~^ hours around these minima 
as illustrated by the inserts in Figure 2. 

The above searches resulted in a collection of period candidates as peaks in the Lomb 
periodogram and/or phase dispersion minima. The light curves were folded with each can- 
didate period to compute an average light curve, which was then compared to the individual 
light curves to check for phase preservation. Only for the period of 32.69 hours were all 25 
individual light curves consistent with the average light curve, while for other periods some 
individual light curves were offset from the average light curve by more than 90°. We con- 
clude that the true period is 32.69 hours. To understand the origin of the period aliases, we 
simulated data using the period of 32.69 hours, with the average light curve described below 
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and the time windows for the observations. The simulated data were able to reproduce the 
period aliases, suggesting they can be attributed to the time windows for the observations 
and the shape of the underlying light curve. 

The average light curve for the period of 32.69 hours, computed with 16 hght curves 
observed in the ultra-deep survey, showed a clear eclipse at phase 0.4 and a minor dip at 
phase 0.95 (Figure 3). The average light curve was compared to all individual light curves, 
including the 18 "good" light curves with more than 100 photons and the other seven light 
curves with fewer than 100 photons (Figure 4). Remarkable consistency in the phase and 
the flux level was found between the average and all individual light curves, except for the 
light curve of ObsID 5340 which was significantly lower than the average. Deep eclipses of 
about eight hours were seen in a number of observations, with steep ingress and egress both 
occurring within one time bin (~1 hour). We thus estimate the eclipse duration to be 8 ± 
1 hours, corresponding to an eclipse full angle of 88° ±11°. Residual emission in the eclipse 
was present in some individual light curves and in the average light curve. The apparently 
slow ingress and egress of the eclipse in the average light curve were artifacts due to the 
averaging of 16 light curves, all of which were not aligned exactly. 

The period estimate can be improved by using the hght curve in ObsID 934, which 
extended the observation time line from one year (~ 268 cycles) to five years (~ 1279 
cycles). For the observations in 2004, the times of the seven minima can be determined to 
better than 30° relative to the average light curve, leading to a fractional error in period of 
AP/P ~ (30°/360°)/268 ~ 0.0003. A slight change in the period docs not change the phases 
much for observations in 2004, but leads to much larger phase changes for the observation 
ObsID 934. By fine tuning the period, the fight curve of ObsID 934 was found to align 
up with the average fight curve to better than 30° at P = 32.688 hours (Figure 5). The 
fractional error in the period is ~ (30°/360°)/1279 ~ 6.5 x 10~^, equivalent to an error 
of 0.002 hours or 7.7 seconds. Thus, the period is 32.688 ± 0.002 hours, and throughout 
this paper we use the following ephemeris: = 0.0 for JD 2450814.49927 + n*(1.36200 ± 
0.00008). 

The source spectrum was extracted for each observation yielding more than 300 counts. 
Six spectra were extracted and fitted by absorbed power- law models (Table 1). The model 
parameters were poorly constrained, and were consistent with each other except for the 
spectrum in ObsID 934. The five similar spectra were combined for fitting to better constrain 
the parameters. This led to uh/W^^ = 3.3 ± 0.5 and F = 3.4 ± 0.3 for the absorbed power- 
law model ixl/dof = 1.06/50), and uh/IO^^ = 9.6 ± 2.6 and T^^ = 0.45 ± 0.4 keV for 
the absorbed multi-color disk model {xl/dof = 1.25/50). The count rates were converted 
to unabsorbed fluxes in 0.3-8 keV with the conversion factor 1 count/ksec = 4.86 xlO"^^ 
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erg/sec/cm^ (i.e. 2.7 xlO^^ erg/sec at a distance of 6.8 Mpc) derived from the combined 
power-law fit. The unabsorbed luminosities were in the range of 0.2-2.3 xlO^^ erg/sec with 
an overall average of 1.3 xlO^^ erg/sec (Table 1). If the combined multi-color disk model was 
used, the conversion factor would be 1 count/kscc = 4.33 xlO~^^ erg/sec/cm^ (i.e. 2.4 x 10^''' 
erg/sec), and the average luminosity was 1.2 xlO^* erg/sec. 

A significant soft excess below 0.5 keV was present in the spectrum in ObsID 934. This 
led to a much lower Uh value (still higher than the Galactic value Hh = 1.6 x 10^° cm^^) 
in the absorbed power law model fit. The spectrum can be fitted with an absorbed power 
law plus multi-color disk model, with uh/W^^ = 1-9 ± 0.7, F = 3.0 ± 0.3 and Tj„ = 64 ± 11 
eV {xt/dof = 1.084/19). The spectral fit was slightly improved as compared to the one 
component models, implying possible presence of a ~ 60eV soft component in 2000 that was 
missing from subsequent observations. However, five other sources, randomly distributed 
over the chip in ObsID 934, exhibited similar soft excesses that were absent in subsequent 
observations, suggesting a possible alternative origin attributable to the uncertainties in the 
calibration of the instrumental soft response in the early observations. 

The spectrum is expected to become harder (softer) at the echpse ingress (egress) due 
to additional absorption by the blocking secondary, and the residual emission in the eclipse, 
presumably from the stellar winds, is expected to be very soft. While the small number 
of photons during the eclipse docs not allow a detailed study of the additional absorption 
and stellar winds, we can at least test whether there is additional absorption and/or stellar 
winds thus hardening and/or softening of the spectrum. The phase for each photon from 
all observations was computed for the period of 32.688 hours. The 340 photons with phases 
0.25-0.50 were considered in the eclipse, and the 4490 photons with phases 0.50-1.25 were 
considered outside the echpse. The "raw" spectra for the two groups were plotted in Figure 
6. Compared to the "non-echpse" spectrum, the "echpse" spectrum exhibited a lack of soft 
photons below 2 keV and an excess of hard photons above 4 keV. A Kolmogorov-Smirnov 
test showed the two spectra were the same at a 0.02 level, indicating the hardening of the 
"eclipse" spectrum was true at a 98% confidence level. To reveal spectral changes associated 
with ingress and egress of the eclipse, we calculated the hardness ratio as (H-S) / (H+S) for 
phase bins of 0.04, as plotted in Figure 7 for three sets of soft and hard bands. The softening 
was present at the minor dip (0 = 0.9) in all three band sets, but with less than 1 — a 
significance. Slight hardening was present at ingress {(f) — 0.25) for band sets S(0.1-1)/H(l- 
8) and S(0.1-1.3)/H(1.3-8), however, it was absent at ingress for band set S(0.1-0.7)/H(0.7-8), 
and the reverse softening is absent at egress for all band sets. Indeed, the hardness ratios 
for most phase bins were consistent with the average hardness ratios within 1 — a error bars 
(Figure 7). 
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3. Optical Observations 

M101-X7 was observed with HST ACS/WFC in the F658N {H^; 2440 seconds) fiher 
on 2004-02-10, and the F435W (B; 900 seconds), F555W (V; 720 seconds), and F814W (I; 
720 seconds) filters on 2002-11-16. The images were downloaded from the STSci archive 
and calibrated on the fly with the best available calibration files. The drizzled images with 
geometric distortion corrected and cosmic rays removed were used to compare with the 
X-ray image from ObsID 6115, in which X7 was off-axis by ~2'. X7 was registered onto 
the optical images with the help of a nearby source, CXO J140339. 3+541827, which was 
positionally coincident with an extremely red bright extended object in a sparse stellar field. 
The positional uncertainty of X7 on the optical images comes from the centroiding errors 
of X7 and CXO J140339. 3+541827 and the ACIS plate scale variation, and is about 0'.'3. 
Due to its large off-axis angle, the X-ray emission of X7 was dispersed to within an ellipse 
of 1'.'4 X 1'.'2 as ovcrplotted on the optical images in Figure 8. 

X7 is located in the center of a super-bubble between the spiral arms as shown on the Ha 
image. This super-bubble shows shell structures with a size of ~ 8"x8" (i.e., ~ 270pcx270pc 
at 6.8 Mpc). The shells exhibit non-concentric irregularities, indicating multiple supernova 
explosions at different off-center locations. Several possible counterparts were found within 
the overplotted eUipse in Figure 8, with the brightest labeled from 'A' to 'J'. To detect point 
sources and compute the photometry, the PSF-fitting package DOLPHOT (Dolphin 2000) 
was run on these images, which resolved some of the labeled counterparts into several stars. 
These counterpart candidates have V = 23.5 - 29.0 mag as listed in Table 2. The V-I colors 
range from -0.7 to 3.3 mag, with an average of 0.7 mag. Such colors indicate the presence of 
heavy reddening commonly seen in star forming regions, or stellar spectral types later than 
GO (V-I = 0.7 for GOV, V-1 = 3.4 for M7V). At the distance of MlOl, dwarf stars later than 
GOV cannot be detected, but G/M supergiants can be detected. 

The counterpart candidate Al was closest to the nominal X-ray position, and had V-1 = 
0.2 + 0.11 mag and V = 24.96 + 0.05 mag, suggestive of an B0-O9V star (My = -4.0 ~ -4.5 
mag) or brighter at 6.8 Mpc (/x = 29.2 mag) in MlOl given the absolute magnitude. To 
search for the spectral types with consistent extinction and reddening, we first calculated 
the extinction Ay for each spectral spectral type above BOV, then the reddening E(B-V) 
= Ay/Rv — Av/3.1 and E(V-I) assuming the standard extinction law (Cardelli, Clayton & 

Mathis, 1989). The search pointed to spectral types of 05-03 (My = -5.7 6.0 mag, V-I 

= -0.47 mag). For such stars, the extinction Ay = 1.5 - 1.8 mag. leading to the reddening 
E(B-V) = 0.48 - 0.58 mag and E(V-I) = 0.57 - 0.69 mag, consistent with the observed E(V-I) 
= 0.67 + 0.11 mag. This reddening corresponded to uh = 2.9 — 3.5 x 10^^ cm^^ (Bohlin, 
Savage, & Drake, 1978), consistent with the uh value derived from the power-law fit to the 
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combined X-ray spectra. 

As a side note, M101-X7 has drawn much attention due to its positional coincidence with 
a supernova remnant candidate MF83 (Matonick & Fesen 1997). Early ROSAT observations 
did not have enough photons and the spatial resolution to determine whether X7 was a point 
source or a diffuse source. If X7 were diffuse emission from shock heated gas in the remnant, 
its high luminosity would require an explosion energy two orders of magnitude higher than 
those of normal supernova and make it a hypernova (Wang 1999). Subsequent ground-based 
observations showed that MF83 is a star formation region with several Hll regions and a 
large ionized shell that is a super-bubble based on its physical properties (Lai et al. 2001). 
If the X-ray emission were indeed diffuse, this shell would be an X-ray bright super-bubble 
whose thermal energy would require either a hypernova or 10-100 supernova in the past 
10^ years. The recent Chandra observations, however, proved that X7 is consistent with a 
point-like source in shape. Furthermore, the presence of periodicity and eclipse in the X-ray 
light curves provides unmistaken evidence for the point-like nature of X7. 



4. Discussion 

MlOl X7 (CXO J140336.0+541925) was observed in 25 Chandra observations span- 
ning five years, and exhibited significant modulations in twelve observations, with apparent 
minima present in seven observations. The Lomb periodogram and the phase dispersion 
minimization were applied in a search for periodicities, revealing a period of 32.688 ± 0.002 
hours for this source. The average light curve folded with this period shows a deep eclipse 
with an eclipse full angle of 9^ ~ 88° ±11°, suggesting that this source is an eclipsing X-ray 
binary in a nearly edge-on configuration. At a distance of 6.8 Mpc, this is the first eclipsing 
X-ray binary discovered outside the Local Group. This discovery can only have been made 
with the superb sensitivity and spatial resolution of the Chandra X-ray Observatory, and 
with the MlOl ultra-deep sequence of observations samphng different orbital phases. 

Residual emission was present in the eclipse, indicating that the X-rays are scattered 
or otherwise processed by, e.g., the strong wind as in the high mass X-ray binaries Cen 
X-3 (e.g., Schreier et al. 1972) and Vela X-1 (Becker et al. 1978). The presence of strong 
winds may lead to a slight over-estimate of the eclipse full angle, and its possible variability 
may have contributed to the mis-alignment of the ingress and egress in some observations. 
The residual emission from stellar winds will make the spectrum in the eclipse appear softer, 
while the additional absorption by the blocking secondary will make it harder. A preliminary 
study, limited by the small number of photons, showed that the hardness of the spectrum 
did not change significantly with the orbital phases, except for slight hardening present in 
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the eclipse. 

Much information can be inferred for this binary from the hght curves. The edge-on 
viewing geometry as suggested by the presence of echpse in the hght curve renders inap- 
phcable the beaming mechanism which enhances the radiation along the polar axis. Thus, 
its average luminosity of Lx=1.3xl0^^ erg/sec constrains the accretor mass through the 
Eddington limit to be larger than 1 Mq. Despite the differences in the average fluxes in 
individual observations, 24 out of the 25 light curves were roughly consistent with the aver- 
age light curve, and only one (i.e., 4%) was significantly lower in fiux than the average light 
curve. This is distinctly different from Be X-ray binaries with typical periods of 100-300 
days, which are in high flux states only briefly when the accretors are near the perigee of 
their highly elliptical orbits. The relatively shorter period, the wide eclipse, and the nearly 
uniform distribution of the high flux states suggest that this binary has a much more compact 
and highly circularized orbit. 

The presence of the 32.69 hour period and the X-ray eclipse place strong constraints 
on the binary properties if we assume the period is the orbital period and the donor fills its 
Roche lobe. The period P can be computed with Kepler's law P — I-ko?!"^ j ^jG{Mx M2), 
in which the separation a can be related to the Roche lobe radius i?cr by R^r — a ■ f{q) — 
a ■ 0.49g2/^/[0.6g2/^ + ln{l + g^/^)] (Eggleton 1983). Here Mi is the accretor mass, M2 is the 
donor mass, and q = M2/M1 is the mass ratio. If the donor is filling its Roche lobe, i.e., 
R2 = Rcr, the period will be determined solely by the accretor mass and the donor mass 
and radius that can be determined empirically by the donor's spectral type. A discrete set 
of donors of the spectral types as tabulated in the Astrophysical Quantities (Allen 2000) are 
tried to obtain this period for an accretor mass range 1-10^ Mq. This period can only be 
obtained for some dwarf stars; all giants and supergiants give much longer periods. 

The donor star can have evolved off the main sequence even if it is not yet a gi- 
ant /super giant star. To consider such donors, we flrst considered main sequence dwarf 
donors spanning a continuous mass range of 0.5-120 Mq, with radii R2 — ro interpolated 
from the empirical mass-radius relation based on dwarf stars of all spectral types. We then 
introduced a helium core mass Mc for each star, using a standard relationship between core 
mass and stellar radius: R2 = ro + 3700M^/(1 -|- M^ + 1.75M^). In the donor mass-accretor 
mass phase plane the viable binary systems with donors of Mc = are shown as the thickest 
line, while the viable binary systems with donors of Mc > populate the regions below and 
to the left of this line (Figure 9). The maximum eclipse full angles Be are computed for the 
viable binary systems with sin{9E/2) = R2/a assuming the inclination angle i = 90°, and 
the observed large 9e (= OEsini) favors those systems with massive donors (Figure 9). Such 
systems have a total mass < IO^Mq, and the radii of the orbits are less than 1 light minute. 
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rendering barycenter corrections unnecessary. 

The accretion rates are calculated for the viable binary systems assuming donors over- 
flow the Roche lobe under conservative mass transfer. Invoking conservation of total angular 
momentum, the accretion rate can be expressed as M2/M2 — N/D. Here the denominator 
can be expressed as 4[1 — (Mc/McY] + 2(1 — g) — (1 + q)dlnf / dlnq for conservative mass 

transfer, and the numerator N can be approximated as "^'^^^^ '^ [ '"^^m^^^ when magnetic 
braking dominates the orbital angular momentum loss. Calculations show that the accretion 
rate decreases for smaller mass ratio g's and larger Mc's; the contours are plotted in Figure 
9 for M2 = 1, 0.1, 0.01, 10^^ 10"*^ and 10"^ in unit of Ix/c^ For a region with g > 1 to the 
right of M2 = Lxjf? 1 D < 0, indicating the breakdown of the adopted expression and very 
large mass transfer rates that likely lead to common envelopes unless the accreted materials 
can be dissipated effectively. 

Similar orbital periods and eclipses have been detected in two well studied X-ray binaries, 
i.e., the intermediate-mass X-ray binary Her X-1 and the high-mass X-ray binary LMC X-4 
as shown in Figure 9. Her X-1 is an eclipsing binary pulsar with a pulse period of 1.24 
seconds and an orbital period of 1.7 days discovered by Uhuru (Tananbaum et al. 1972). 
Its companion, HZ Her, is a ~ 2Mq A/F star (Bahcall & Bahcall 1972) fiUing its Roche 
lobe in a compact circular orbit. The X-ray light curve for Her X-1 consists of alternating 
main on-state, off-states, and short on-states with a super-orbital period of 35 days, which 
is believed to be caused by the precession of a tilted accretion disk, and has been studied by 
almost all X-ray missions (cf, Shakura et al. 1998). LMC X-4 was first discovered by Uhuru 
(Giacconi et al. 1972), and was identified with an 07 III-V star exhibiting a binary period of 
1.4 days (Chevaher & Ilovaisky 1977). The X-ray pulsations were discovered with a period 
of 13.5 seconds, the Doppler variations of which revealed the system as a neutron star in 
circular orbit with a 17 Mq donor nearly filling its Roche lobe (Kelley et al. 1983). The X- 
ray light curve for LMC X-4 exhibits echpses, occasional flares, and a 30.5 day super-orbital 
modulation as prototyped by Her X-1 (Lang et al. 1981). 

MlOl X7 may be an MlOl analog for Her X-1 or LMC X-4 given the similarities in the 
orbital period and the presence of echpses. There is, however, no evidence for a super-orbital 
period of about 30 days for MlOl X7, and the pulse period cannot be detected given the 
low count rate of about 0.005 cnt/sec. Furthermore, its duty cycle for the low flux state is 
~ 4%, while the duty cycle for the low flux state (off-state) is >50% in Her X-1 and LMC 
X-4. In addition, its eclipse full angle is much larger than those for Her X-1 and LMC X-4 
(both ~ 50°), suggesting a higher inclination angle and/or a more massive donor than for 
Her X-1 or even LMC X-4 (Figure 9). Note, however, that the uncertainties in the duty cycle 
and eclipse angle are large, given the small number of photons collected and the temporal 
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sampling. This makes direct comparison with Her X-1 and LMC X-4 difficult. 

The donor and accretor masses can be further constrained by monitoring and identifying 
the donor in the optical. Our analyses of the HST ACS/WFC archive observations have 
resulted in several candidates for the donor, which span a large range of magnitude and 
mass. The donor mass, thus the primary mass, can be further constrained if the donor can 

be identified with a single counterpart. For example, were the 05-03 star Al (60-120 Mq) 
identified as the secondary, the binary system would be constrained to the lower right corner 
in Figure 9, with the primary mass lower than 10 Mq. This identification can be achieved by 
photometric monitoring of these counterparts over an orbital period, because the optical light 
curve should exhibit eclipse and ellipsoidal modulations given the edge-on viewing geometry 
of the binary. These observations can be carried out with either HST or large ground-based 
telescopes with adaptive optics such as Gemini-North in the north hemisphere. 
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Fig. 1. — Scarcli for periodicities for M101-X7 with the Lomb-Scargle method. The peri- 
odogram was computed for the 18 "good" light curves. The power 25.8 corresponds to the 
probabihty of 10~^° for the periods to originate from random fluctuations in the absence of 
true signals. 
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Fig. 2. — Search for periodicities for M101-X7 with the phase dispersion minimization 
method. The phase dispersion was calculated for the 18 "good" light curves, first with 
trial period steps of 0.045 hours, then with finer steps of 2 x 10^"^ hours around the local 
minima. The three inserts show the finer searches around 22.82 hours, 32.72 hours, and 
44.87 hours as examples. 
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Fig. 3. — The light curves folded with the period of 32.69 hours for the 17 observations in the 
MlOl ultra-deep Chandra survey. The thickest histogram is the average light curve from 16 
light curves excluding ObsID 5340, which is significantly lower than the average light curve. 
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Fig. 4. — The individual light curves from the MlOl ultra-deep survey in comparison to 
average light curve for the period of 32.69 hours. The light curves with a "*" suffixed in 
key are the 18 "good" light curves used to search for the period. 
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Fig. 5. — Tlie liglit curve for ObsID 934 in comparison with the average hght curves for 
shghtly different periods. The hght curve are mostly consistent with the average hght curve 
for the period of 32.688 hours (thick sohd) by ahgning the echpse ingress and the minor dip. 
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Fig. 6. — The energy distributions for photons within the echpses and outside the echpses 
from all 25 observations. The coadded spectrum within the eclipses exhibited a lack of soft 
photons below 2 keV and an excess of hard photons above 4 keV as compared to the coadded 
spectrum outside the eclipses. 
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Fig. 7. — The hardness ratios as a function of phases. The hardness ratio was calculated 
as (H-S)/(H+S). Phase bins of 0.04 were used to ensure at least 40 photons per bin. The 
straight lines are the average hardness ratios for three band sets and their errors. 
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Fig. 8.— (a) The HST ACS/WFC image for XT. The 9" x 9" box encloses the irregular 
shell structure of the super-bubble. The ellipse centered at the corrected X-ray position 
of X7 encloses 95% of the source photons as reported by WAVDETECT. (b) The HST 
ACS/WFC B image for X7. The ellipse is the same as in the Ha image. The letters label 
the concentrations of bright stars, (c) The HST ACS /WFC V image for X7. (d) The HST 
ACS/WFC I image for X7. The overplotted ellipse has a size of 1'.'4 x 1'.'2, while the corrected 
X-ray position of X7 has an error of about O'.'S that only encloses stars A1-A3. 
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Fig. 9. — Possible binary systems with the 32.69 hour period in the donor mass-accretor 
mass phase plane. The thickest line stands for the viable binaries with donors interpolated 
from the dwarf stars of all spectral types. The region below and to the left of this line 
are populated by viable binaries with different He core mass, and the dotted lines denote 
Mc=0.07 and 0.13 M©. The thin solid lines denote the accretion rates M2=l, 0.1, 0.01, 10"^ 
10~^ and 10~^ Lxjf?-, respectively. The systems to the right of M2 = Lx/c^ may develop 
common envelopes unless the accreted materials can be dissipated effectively. The dashed 
lines denote the maximum eclipse full angles 9e = 80°, 60°, 40° and 20°, respectively. 
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Table 1. Individual observations for CXO J140336.0+541925^ 
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^The columns are (1) ObsID, (2) MJD, (3) exposure time in seconds, (4) Off-Axis Angle of 
the source in arcseconds, (5) wavdetect detection significance, (6) net counts derived from aperture 
photometry, (7) average net count rate in ct/ksec, (8) Lx(0.3 — 8/cey) in 10^^ erg/sec, (9) absorbing 
column density in unit of 10^-*^ cm~^ in the power law model, (10) the power law photon index, (11) 
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the reduced and degrees of freedom, and (12) the null hypothesis probability. 

''Observations showing apparent variations are suffixed with a an extra '*' is suffixed if the 
observation contains an apparent minimum; a third '*' is suffixed if steep ingress and egress of a 
deep eclipse are present. 

^1 ct/ks is converted to 2.7 x 10^^ erg/sec in 0.3-8.0 keV using the combined power-law fit, and 
a distance of 6.8 Mpc. The conversion factor is 1 ct/ks = 2.4 x 10^^ erg/sec if using the combined 
multi-color disk fit. 

'^Thc absorbed multi-color disk model is fitted in 0.3-8.0 keV, with Ufj specified to be no less 
than the Galactic value of 1.6 x 10^° cm~^. The second parameter ('T") is the disk inner edge 

temperature Tin- 

^Thc power-law model is fitted to the band 0.7-8 keV by freezing Ufj = 3.3 x lO^"' cm~^ derived 
from the combined fit, to exclude the spurious soft excess below 0.5 keV. 



-25- 



Table 2. Magnitudes for X7 counterpart candidates 
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G 
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^The listed magnitude is for the counterpart D, which was not resolved on the image. 
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